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The exact and analytical low energy spectrum of the ferromagnetic Kondo lattice model on 
a complete graph extended with on-site repulsion U is obtained in terms of a dynamical spl(2,l) 
supersymmetry in the limit of infinitely strong Hund's coupling (double exchange). Furthermore, we 
show that for the particular value of (7 = Jh /2 the supersymmetry is not constrained to infinite Jh /t 
only and we calculate the energy including the /Jh corrections analytically and we give numerical 
evidence which suggest that the Kondo Hamiltonian is itself supersymmetric for any Jh /t. On a N 
site graph, the ferromagnetic ground state is realized for 1 and N+1 electrons only. In the leading 
order in the value of the core spin, the quantum and semiclassical spectra are identical. 
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The recent experimental activity on manganites, 
Ij&i-xAx MnOs (where A=Ca,Sr or Ba), a family in- 
tensively studied due to its colossal magnetoresistance 
[1] , has stimulated the interest of theorists for this com- 
pound. While the orbital degrees of freedom certainly 
cannot be neglected for the real material [2], it is be- 
lieved that some aspects of the physics of the transition 
metal oxides can be revealed by considering the Kondo- 
lattice Hamiltonian TYkl = T + TL-mt , where the kinetic 
part 
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Uj > 0, 



(1) 



describes the hopping of electron on a lattice {a =1,1 is 
the spin index), and the interaction part reads 
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(2) 



Here Jh stands for the intratomic ferromagnetic ex- 
change (Hunds coupling) between the conduction elec- 
trons and localized core electrons [3], <t„^ denotes the 
vector of Pauli matrices, and Scj is the spin operator of 
the localized core electrons with spin Sc- We also in- 
clude the on-site Coulomb repulsion U > Q between the 
electrons, with U of the order of Jh- While here we are 
interested in Jh > 0, the Kondo-lattice Hamiltonian has 
been extensively studied for Jh < 0, as describing heavy 
fermion systems [4]. 

Because the Hund coupling is typically larger than the 
hopping (in manganites Jh « 1 — 2eV and t 0.1 — 0.5 
eV), the energetically unfavorable low spin states are 
usually neglected and we get the quantum double ex- 
change, with a rather complicated Hamiltonian [5,6]. 
The next usually used approximation neglects the quan- 
tum spin fluctuation of the Sc core spin replacing it 
by classical variables (spherical angles 6 and (j)). These 
two approximations lead to the double-exchange model 
[6,7] with Hamiltonian H^c = ^^Ujfjfj : which de- 
scribes noninteracting spinlcss fcrmions (charges) mov- 
ing in a disordered background of classical spins ttj = 



Uj [cos(e',/2)cos(6lj/2) +sin(e',/2)sin(6'j/2)e'('^'-'^^)] . 
The charges can freely propagate provided the core spins 
are aligned and therefore ferromagnetism is favored. The 
main effect of finite Jh is to introduce antiferromagnetic 
exchange between the core spins, which will hinder the 
free propagation of the holes, resulting in a competition 
between ferromagnetic and antiferromagnetic ordering. 
To check the scenario presented above, different numeri- 
cal methods, like exact diagonalization, quantum Monte 
Carlo, density matrix renormalisation group and dynam- 
ical mean field are applied [8-11]. 

Here we will consider the model where the hopping is 
infinitely long ranged, i.e. complete graph with tij = 
t{l — Sij). The numerical diagonalizations show that 
the spectrum separate into well defined bands for large 
enough Jh- A typical example of the lowest band is 
shown in Fig. 1 in the case of a small 4-site cluster where 
states at the right upper corner belong to the next band. 
It appears that with long ranged hopping, the lowest 
band is surprisingly simple in the limit Jn/t — > cxd, and 
U independent with huge degeneracies. The distribution 
of the eigenvalues is analogous to the spectrum of the 
U/t — > oo Hubbard model (t-model) on a complete graph 
[12,13] and, as shown later, this peculiar spectrum re- 
flects the same underlying spl(2,l) dynamical supersym- 
metry [12] (supersymmetry in the sense of odd and even 
number of fcrmions mixed in the same irreducible rep- 
resentation). For finite Jn/t the degeneracies are lifted 
except in the special case U = J///2 where some states 
with the same total spin value S remain degenerated as 
it can be realized in the qualitative comparison of Figs la 
and lb. To understand this behavior: (a) we introduce 
Schwinger bosons to describe the spin degrees of freedom, 
and the form of the effective strong coupling Hamiltonian 
becomes remarkably simple, especially for U = Jh/'^', (b) 
the effective strong coupling Hamiltonian is then solved 
exactly using dynamical supersymmetry and the spl(2,l) 
graded algebra; (c) we show that in the limit Sc ~^ oo the 
semiclassical and quantum model give the same energy 
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FIG. 1. The low energy spectrum of the ferromagnetic 
Kondo lattice on a 4 site complete graph with Sc = 1/2 and 
iVe = 2 for [/ = Jh/2 (a) and [/ = (b). Each point is highly 
degenerate. The solid straight lines show the energy of the 
effective Hamiltonian as given by Eqs. (9)-(12). 



spectrum, and that the scmiclassical limit can be pc;r- 
formed already at the level of the Kondo lattice Hamil- 
tonian. 

a. Effective strong coupling Hamiltonian. In the 
atomic limit (t = 0), the different lattice sites decouple. 
An empty site has energy 0; one electron can form with 
the core spin either the high {Sc + 1/2) or low {Sc — 1/2) 
spin state, with E = -JhSc/2 and E = Jh{Sc + l)/2, 
respectively; finally two electrons results in a state with 
energy E = U and spin Sc- This can be summarized by 
representing the different states using auxiliary fermions 
/ and d and the spins by Schwinger bosons b„ [14]: 
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y/2Sc 



- ^ (b^ ft 



where and S are the creation operators of high and 
low spin states respectively while the spin operators in 
terms of Schwinger bosons are S^ = (nj^ — nJ|)/2, 

S'^ = and SJ ~ ^]j,^jT' "^^^ anticommutation 

relation for electrons require the following constraint to 
be satisfied at each site 



E 



n]c, - n^j + n'^j = 2Sc 



In this representation the interaction part becomes diag- 
onal 

^int = ^ E [^^- + - ^-""J - ^^"^'] +^ E "f^J- 

j j 

Choosing U = Jh/"^ we can eliminate the four fermion 
term njUj. The effective strong coupling Hamiltonian 
is then the expansion in t/Jn around the atomic limit 
and is obtained using a canonical transformation. As a 
first step, following Ref. [15], we decompose the kinetic 
Hamiltonian into a part conserving the number of d and 
/ fermions (To), and a part where the final and initial 
state differs by a large energy Jh{Sc + 1/2) 

^^ = -E^(M.T-M.i)4/.' 



7j. The effective Hamiltonian is 



so that T = To+Tj 
then given by HeS = Hmt+To 

0{t^/J'^), and in the lowest energy subspace, where we 
keep the / fermions only, it reads 



/[(5e + l/2)JH] + 



eff 



Ju^eSc 



2Sc 



jH{2Sc+iy 



The first order term proportional to t in the expansion 
is equivalent to quantum double-exchange Hamiltonian 
derived by Kubo and Ohata [5] and by Miiller-Hartmann 
and Dagotto [6]. Let us note here, that the procedure 
can be repeated for general U, and it would give us a 
more complicated t^ / Jh term [16]. 

h. The model on a complete graph. In this case, we 
can write the Hamiltonian (3) as 

Weff = NcS. - P^cPa (4) 



2Sc + 



2t' 



Jh{2Sc + If 



1)3 E [(a^^c + iVe/2) 6aP - SfTc^ff] F^, 



a,fS 



where = t — JhSc/2, and we introduced the fcrmionic 

operators = J2j bjcfj and = J2j bjafj- with the 
nonvanishing anticommutation relations 

{Fl,F0}=YSc^0 + Sac,0, (5) 

where the spin operators S = J2j satisfy the usual 

[S+,S-]=2S', [S',S^] = ±S^, (6) 
su(2) spin algebra, and with the Y operator defined as 

Y = {Sc + l)N-Nc/2, [S,y]=0. (7) 
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Finally, the commutation relations between the fermionic 
and bosonic operators read: 

[F„,S] = i^<T„0F^, [F^,Y] = -FJ2, (8) 

with their conjugate. The set of relations (5)-(8) define 
a spl(2,l) graded algebra [17,18]. 

The rank of spl(2,l) is 2, and we can define two linearly 
independent Casimir operators [17]: 

C, = S^-Y^-\Y.{FiF^-F^Fl), 

a. 

\ / a/3 

In general case the irreducible representation \Y, S] is 85* 
dimensional and contains the {S,Y), (5 — 1/2, F + 1/2), 
(5- 1/2, y- 1/2) and {S-1,Y) spin multiplets. Special 
cases concern the irreducible representation \Y = S, S] 
which contains only the {S,Y) and {S - 1/2,Y + 1/2) 
spin multiplets (dimension 45+1) and the irreducible 



representation [Y. S ~ 1 /2] which do not contain the 
(S — l,y) spin multiplet. Applying operators F^ and 
F\ we can walk between the different («?, F) spin mul- 
tiplets within an irreducible representation. The eigen- 
values of the Casimir operators Ci and C-i are S"^ — Y"^ 
and (5^ — y^)y in the irrep \Y, S], respectively, while the 
eigenvalues of operators and Y are S'{S' + 1) and Y' 
for an (S",y') multiplet. 

We can find the eigenvalues of the effective Hamilto- 
nian using the representations of the spl(2,l) superalge- 
bra. The Hamiltonian conserves both the electron num- 
ber and total spin: \H, Y] = [H, S] = 0, however it does 
not commute with the F's. On the other hand, since 
EaF^Fa = S^-r'+Y-C2, the (5-1/2, y-1/2), {S,Y), 
(5-l,y), and (S -1/2,Y +1/2) multiplets of the irre- 
ducible representation [Y, S] are eigcnstates of F^^F^ 
with the eigenvalues 2y— 1, Y+S, Y — S, and 0, respec- 
tively. Similarly, the Cz appears in the /Jh corrections 
in the Hamiltonian. Finally, the effective Hamiltonian 
turns out to be a combination of C2, C3, and Y, and 
the states \{S,Y)[y',s']} are eigenstates with (increasing) 
energies: 



(9J.2 omax 1 1 \ / cmax \ A4-2 q( q \ T\ +3 



(9J.2 cmax _ C \ / cmax _ c\ +3 



(10) 



V, - ... - (. . £9^) (« - mi) . o,il,, 

E[Y-i/2,s+i/2]{S,Y) = N,e_, (12) 



where 5""'^'^ = NSe + Ne/2 = {2Sc -M) - y is the max- 
imum total spin. In Fig. 1 these energies are compared 
with those of the Kondo lattice for U = Jh/2 at left. It 
appear that the computed correction is the correct one 
but for t/Jji > 0.05 higher terms are necessary. It is in- 
teresting to note that the E[y-i/2.s+i/2](S,Y) = S-Ne 
is the eigenvalue for any Jh- To obtain this solution it is 
essential that the effective Hamiltonian can be expressed 
using the operators of the spl(2,l) superalgebra which, 
for finite Ju/t, is possible for U = Jh/2 only. 

Let us now address the question of the degeneracy of 
each level. For one site, the states with S = S,- and 
5c -1-1/2 form the irreducible representation [5c -1-1/2, 5c -|- 
1/2]. With N sites, the de generacy M^y is determined 
by the number of times the irrep [Y, 5] is contained in 
[5c-t-l/2, 5c-hl/2]^. This is determined from the branch- 
ing rule [17], leading to the following recursion relation: 

S-\S-y-Sa\ 

-'"[r,S-|-l/2] — ^^^[Y+y-S^-l,S-s+l/2Y 
y={0,l/2} s=-Sc-y 

with M^[YlN(Sa+i/2),s] = h;N(Sa+i/2) as boundary con- 
dition. This recursion formula along with the energy def- 



inition in Eqs. (9)-(12) allows an iterative procedure to 
build the energy density distribution of the model. In 
Tab. I we show as an example the spectrum of the 4 site 
cluster. 

For Ne = l the (5, Y) spin multiplet of the [Y+1/2, S+ 
1/2] is missing and the groimd state is the highest spin 
state in the [Y,S] irrep. For 1 < A^g < the I/Jh cor- 
rection makes the lowest energy state to be the singlet 
in the [Y + 1/2, 1/2] irrep, and the low energy spectrum 
behaves as 5(5 -|- 1), like in the antiferromagnetic infinite 
range Heisenberg model. For Ng > N we can use that 
E{Ne, t, Jh) = E{2N - N,,, -t, Jh) + U{Ne - N), and 
from this we obtain that when A^g = A' -h 1 the ground 
state is again the highest spin state (a similar situation 
is encountered in the Hubbard model [19]). Finally, for 
A^ -|- 1 < A^e < 2N the ground state is highly degener- 
ate and the corresponding energy can be obtained from 
Eq. (12) with the appropriate change. To summarize, the 
highest spin multiplet is the non degenerate ground state 
for Ne = 1 and Ng = N + 1 only, cases referred to true 
ferromagnetism. 

c. Semiclassical limit. In the limit 5c ^ oo we can 
replace Sj by 5c (sin Oj cos , sin 9j cos 4>j , cos 9j ) in the 
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Kondo Hamiltonian (2) with J = Jh{Sc + 1/2) fixed, 
and we end up with an noninteracting model [20]. The 
one-particle spectrum comes out as follows: two N — 2 
degenerate states with energy e± = t ± J/2 and four 
nondegenerate states with energies 

esrs. = + '-l^j2+NH^-2s2Jt^^ , 

where si and S2 are ±1. The lowest energy state corre- 
sponds to Si = — 1 and S2 = —1, and for Jh ^ t the 
energy is linerly decreasing with 5, i.e. the fermions can 
freely propagate when the spins are parallel. In the next 

state (e ^) the tendency is reversed: energy is higher 

for larger 5*. The absence of the true ferromagnetic state 
can be traced to the cancelation of the contributions lin- 
ear in S for Np, > 2. This does not happen in a model 
with less pathological one-particle spectrum, e.g on the 
cubic lattice the ferromagnetism is realized for a wide 
range of hole concentration [8]. Filling the one-particle 
levels, the lowest part of the spectrum for 1 < iVg < 

electron is e -I- e |_ -I- {Ne — 2)e_, e__ -I- {Ne — l)e-, 

e \- + {Ne—l)e-, and Ne£-. These energies are equal, up 

to corrections 0{t^ / J"^, ^/Sc), to the energies (9)-(12) of 
the {S, Y) spin multiplet in the [Y + 1/2,8+ 1/2], [Y, S], 
[Y,S + 1], and [Y - 1/2, S + 1/2] irreps, respectively. 
This way we can establish a one-to one correspondence 
between the semiclassical and quantum spectra for this 
particular model. 

To conclude, we have shown that the effective strong 
coupling limit of the Kondo lattice model on a complete 
graph can be described by the spl(2,l) graded algebra 
and exhibits a dynamical supersymmetry. In contrast to 
Hubbard model, the supersymmetry is not limited to the 
t-model only, but holds also in the next order {t^ / Jh) for 
a special value of on-site repulsion. This, and the numeri- 
cal diagonalization on small clusters suggest that not only 
the strong coupling limit, but the Kondo Hamiltonian it- 
self is supersymmetric. We also show that for — > oo 
the spectrum of quantum model can be derived from one- 
particle states of the semiclassical Kondo lattice model. 
Finally, the ferromagnetic ground state is reduced to the 
case when we have only 1 or N-|-l electron, in contrary 
to the generally accepted expectations. 
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TABLE I. The spectrum and multiplicities of the quan- 
tum double exchange (Jn/t — > -|-oo) for 5c = 1/2 and N = A. 
Here E' = E — NeS-. As an example, we have underlined 
the four spin multiplets of the [Y = 5,S = 2] irreducible 
representation, which has multiplicity 8. 
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